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In this paper, we address the design of high spectral-efficiency Barnes-Wall (BW) lattice codes 

which are amenable to low-complexity decoding in additive white Gaussian noise (AWGN) channels. 

c/3 ' We propose a new method of constructing complex BW lattice codes from linear codes over polynomial 

rings, and show that the proposed construction provides an explicit method of bit-labelling complex BW 
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upper bound on its error performance. We show that the SBWD is powerful in making correct decisions 
well beyond the packing radius. Subsequently, we use the SBWD to decode lattice codes through a 
novel noise-trimming technique. This is the first work that showcases the error performance of SBWD 
in decoding BW lattice codes of large block lengths. 

Index Terms 

Barnes-Wall lattices, lattice codes, low-complexity lattice decoders. 

*The authors are with the Department of Electrical and Computer Systems Engineering, Monash University, Melbourne, 
Australia-3168. Email:harshan.jagadeesh@monash.edu, emanuele.viterbo@monash.edu. ^The author is with the Department 
of Communications and Electronics, Telecom ParisTech, Paris, France, Email:belfiore@enst.fr. Parts of this work are in 
the proceedings of IEEE International Symposium on Information Theory (ISIT) 2012 held at Cambridge, MA, USA, and 
International Symposium on Mathematical Theory of Networks and Systems (MTNS) 2012 held at Melbourne, Australia. 

January 9, 2013 DRAFT 



I. Introduction 

Ever since random coding schemes were demonstrated to approach the capacity of additive 
white Gaussian noise (AWGN) channels [H], enormous research has taken place to find structured 
coding schemes which can accomplish the same job. The need for structured coding schemes is 
to facilitate simpler analysis of the code structure and to achieve reduced complexity in encoding 
and decoding. A well known method of obtaining structured codes is to carve out a finite set 
of lattice points from dense lattices [|2l-[|5l. Such codes are referred to as lattice codes, and are 
usually obtained as a set of coset representatives of a suitable quotient lattice. Further, the lattice 
codes have the advantage of inheriting most of the code properties from the parent lattice, and 
as a result, the choice of the lattice is crucial to the performance of the code. 

A. Motivation and contributions 

In this paper, we are interested in carving lattice codes from Barnes-Wall (BW) lattices [lH, 0. 
Our goal is to construct efficient BW lattice codes of large block lengths which work with low- 
complexity encoders and decoders. In particular, efficient low-complexity decoders for complex 
BW lattices are readily available in [[T3l , [[T5l . Therefore, if lattice codes from complex BW 
lattices are employed for communication over AWGN channels, then the decoders of [|T3l , [|T5l 
can be used to recover information with low computational complexity. 

In |fT3l , the authors have proposed two low-complexity implementations of the bounded 
distance decoder for BW lattices, namely (?) the sequential bounded distance decoder, and (//) the 
parallel bounded distance decoder. Inspired by the parallel decoder in [[T3l , list decoders for BW 
lattices have been recently proposed in [[TSl . We note that the parallel decoders of [[131 and [[151 
have low-complexity only when implemented on sufficiently large number of parallel processors. 
If the above decoders are implemented on a single processor, then the complexity advantages are 
lost, and specifically, the complexity of the list decoder grows larger than that of the sequential 
decoder in [[T3l . Since we are interested in lattice codes of large block lengths, we focus on the 
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sequential bounded distance decoder which seems more suitable for implementation (see Section 
IV-BI for more details on the complexity advantages of sequential bounded distance decoder over 
the list decoder). The sequential decoder in [[131 was proven to correct any error up to the packing 
radius. However, the possibility of correct decoding is not known when the received vector falls 
outside the bounded decoding ball of packing radius. In a nutshell, the exact error performance 
of the decoder is not known. The existence of this low-complexity decoder has motivated us to 
study its error performance, and use it to decode BW lattice codes. We refer to this decoder as 
the sequential BW lattice decoder (SB WD). The contribution of this paper on the construction 
and decoding of complex BW lattices are given below. 

1) We introduce Construction A' of lattices which enables us to generate some well structured 
lattices from linear codes over polynomial rings [124] . As an immediate application, we 
apply Construction A' to obtain BW lattices of dimension 2™ for any m > 1. The proposed 
method is yet another construction of BW lattices (shown in Section HIH) and shows a new 
connection between codes over polynomial rings and lattices. We show that the proposed 
construction provides an explicit method of obtaining and bit-labelling complex BW lattice 
codes. 

2) We study the error performance of the SBWD in AWGN channels. Since the error per- 
formance of the SBWD depends on the error performance of the underlying soft-input 
Reed-MuUer (RM) decoders, we study the error performance of the soft-input RM decoder 
as used in the SBWD. First, we use the Jacobi-Theta functions [^261 to characterize the 
virtual binary channels that arise in the decoding process. Subsequently, we study the noise 
statistics in the algorithm, and provide an upper bound on the error performance of the 
soft-input RM decoders. Through computer simulations, we obtain the error performance 
of the SBWD, and show that the decoder is powerful in making correct decisions well 
beyond the packing radius [|25l (see Table I in Section |V] for the effective radius of the 
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SBWD decoder). This is the first work that showcases the error performance of SB WD in 
decoding BW lattices of large block lengths. 
3) To decode the lattice code in AWGN channels, we employ the SBWD along with a noise 
trimming technique, wherein the components of the received vector are appropriately scaled 
before passing it to the SBWD. With the noise-trimming technique, the SBWD is forced 
to decode to a codeword in the code which in turn improves the error performance. We 
refer to this decoder as the BW lattice code decoder (BWCD). We obtain the bit error rate 
(BER) of the BWCD for codes in complex dimension 4, 16, and 64, and show that the 
BWCD outperforms the SBWD by 0.5 dB. 

B. Prior work on Barnes-Wall lattices 

The BW lattices [6J is a special family of A^-dimensional lattices that exist when A^ is a 
power of 2. These lattices were originally discovered as a solution to finding extreme quadratic 
forms in 1959 [|6l. Only in 1983, the now well known connection between BW lattices and 
Reed-MuUer codes was discovered by dHl. This connection is found in several works [|9l, [[TOl , 
lfT4]| in different forms. Other than its construction through Reed MuUer codes, the generator 
matrices of the BW lattices are also known to be obtained though Kronecker products [jTll . [[T3l . 

In 1989, G.D. Forney proposed a low-complexity bounded distance decoding algorithm for 
Leech lattices [fTSl . As a generalization, in the same paper, a similar algorithm has been shown 
to be applicable in decoding all Construction D lattices. As BW lattices are known to be 
obtained through Construction D [Q, bounded distance decoders for BW lattices were known 
in principle since lITSl . Only in the 1990's, explicit bounded distance decoders for BW lattices 
were implemented for dimension up to 32, and numerical results on the error performance were 
reported [[T9l , [|20l . In 2008, Micciancio and Nicolosi liT3l have proposed two low-complexity 
implementations of bounded distance decoder for BW lattices, namely (/) the sequential bounded 
distance decoder, and {ii) the parallel bounded distance decoder. If A^ = 2™ denotes the dimension 
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of complex BW lattice, the worst-case complexity of the decoders has been shown to be 
O (Mog^(A^)) and O (log^(A^)) for the fully sequential decoder and the fully parallel decoder, 
respectively. For the fully sequential decoder, the algorithm is assumed to be implemented 
on a single processor, whereas for the fully parallel decoder, the algorithm is assumed to be 
implemented on N"^ parallel processors. Inspired by the fully parallel implementation in [IT3l|, 
list decoders for BW lattices have been recently proposed in IfTSl where the list decoder outputs 
a list of BW lattice points within any given radius from the target vector. The complexity of the 
list decoder is shown to be polynomial in the dimension of the lattice, and polynomial in the list 
size, which is a function of the Euclidean radius. Note that the SBWD exploits the Construction 
D structure of BW latices as a multilevel code of nested binary Reed-MuUer (RM) codes, and 
decodes each RM code through a successive interference cancellation technique. On the other 
hand, the list decoder does not exploit construction D structure of BW lattices, and hence, does 
not need the support of any soft-input RM decoders. 

The rest of this paper is organized as follows: In Section III we provide a background on 
lattice constructions from linear codes. In Section Illl we introduce Construction A' of complex 
BW lattices. In Section HVl we study the error performance of the SBWD, while in Section IVl 
and Section |VIl we use the SBWD to decode the BW lattice code. Finally, in Section IVII[ we 
conclude this paper and provide some directions for future work. 

Notations: Throughout the paper, boldface letters and capital boldface letters are used to 
represent vectors and matrices, respectively. For a complex matrix X, the matrices X^, 3fi(X) 
and S(X) denote the transpose, real part and imaginary part of X, respectively. The set of 
integers, real numbers, and complex numbers are denoted by Z, M, and C, respectively. We use 
i to represent y/^. For an n-length vector x, we use Xj to represent the j-th component of x. 
Cardinality of a set S is denoted by \S\. Magnitude of a complex number x is denoted by \x\. 
The number of ways of picking n objects out of m objects is denoted by C^. The symbol [-J 
denotes the nearest integer of a real number, and we set \a + 0.5J = a for any a E 7j. Finally, 
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we use Pr(-) to denote the probability operator. 

II. Background on Lattice Construction using Linear Codes 

A complex lattice A over Z[i] is a discrete subgroup of C" fl9|. Alternatively, A is a TL\i\- 
module generated by a basis set {vi, V2, . . . , v„, | Vj G C"} as A = < YTi=\ 'io^i I ^Q'i ^ ^['^] f • 
It is well known that dense lattices can be obtained via binary linear codes [jOJ. Depending on 
the structure of the underlying linear codes, lattice construction can be categorized into different 
types. In this section, we recall two well known constructions for the case of complex lattices [|9l. 

Construction A: 

Definition 1: A complex lattice A is obtained by Construction A from the binary linear code 
C if A can be represented as 

A = (l + i)Z[^]"©£o, (1) 

where £0 = {^(c) | Vc G C} C A is a lattice code obtained by the component- wise mapping 

t/' : F2 ^ Z[z] given by V'(O) = and ^(1) = 1 on the alphabet of C, where F2 = {0, 1}. 

Construction D: 

Definition 2: A complex lattice A is obtained by Construction D from a family of nested 
binary linear codes Cm_i 3 C.m-2 ^ • • • ^ Ci 3 Co if A can be represented as 

A = (1 + i)'^Z\if © (1 + iT-^C^-i © ■ ■ ■ © (1 + z)£i © £0, (2) 

where Cj = {^(c) | Vc G Cj} is obtained by the component-wise mapping ^ : F2 — )• Z[i] given 
by ?/;(0) = and V^(l) = 1 on the alphabet of Cj. 

A BW lattice can be obtained via construction Z) as a Z[i] lattice as follows [[TOl . Suppose we 
want to construct the complex lattice BW2m of dimension 2"^ where m > 1, let TZAi{r,m) be 
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the binary Reed-MuUer (RM) code (Sec. 3.7, Ch. 3, HH) of length 2"* and of order < r < m. 
Then, BW2"^ can be constructed as 

{m-l "j 

(1 +i)™a+ ^(1 + i)"^(c,) I Vc^ e 7^Al(r,m),Va G Zfz]^'" I (3) 

where ^(■) is as given in Definition [2l For notational convenience, we also write (|3]) as 

m— 1 

BW2^ = (l + ^)™Z[^]2'"©0(l+^)'■7^7W(r,m). (4) 

r=0 

This method generates BW2m as a multi-level structure of nested RM codes and hence it falls 
under Construction D [Q. 

Generalized construction A lHHl lllTll ; 

Apart from Construction D, the BW lattice codes can be obtained by the generalized con- 
struction A. For the complex BW lattice BW2-^, let Gbw ^ C^^^ denote the generator matrix 
in the triangular form, where the rows {g^^, §2, • • • ,gAr} of Ggw forms a basis set of BW2^, 
where N = 2"^. Let di, d2, . . . , d^ represent the diagonal elements, where dj = (1 + z)™^ for 
some integer nij > 0, and d = (1 + ^ynaxm^ p^^. j-j^j^ lattice construction, one can easily map 
binary data to lattice points as follows: 

1) Bit Labelling: Map log2(-Lj) information bits to aj E Z[i]/pjZ[i] where pj = j- and Lj 
is the cardinality of Z[z]/j>jZ[z]. 

2) Encoding: Using {ai, 02, ... , a^}, a lattice point is obtained as J2j=i '^jSj- 

3) Shaping: Since A = dZ[i]^ + C, a lattice point within £ can be obtained as x = x mod 

dZ\t]^. 

Motivation for Construction A': 

In the bit-labelling step above, binary digits have to be mapped to the symbols of Z[z]/pjZ[z]. 
Some of the well-known bit-labelling methods include gray-mapping and set-partitioning based 
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methodslll Unlike the case of real integer lattice, Z[i]/pjZ[i] is an arbitrary subset of Z[i], and 
bit mapping to Z[i]/pjZ[z] is not straightforward unless the set of representatives for Z[z]/j>jZ[z] 
is chosen with good shaping property. Through Construction A', we facilitate bit-labelling on 
complex integers by using the truncated binary expansion of the elements of Z[i]/pjZ[i] over 
the base 1 + 2 [|23l . With this, the bits labelled on Uj E Z[i]/j9jZ[i] are nothing but the bits in the 
truncated binary expansion of Uj. To assist the bit-labelling step, we use polynomial rings over 
F2 to represent the elements of Z[i]/pjZ[i]. For the encoding step, we use a linear code over 
polynomial rings, and obtain the lattice points as embedding of the codewords a of linear code 
into the Euclidean space. Finally, for the shaping step, we propose an appropriate mapping on 
Z[i] which provides a label code with appropriate shaping property, i.e., we explicitly provide a 
method of bit-labelling complex BW lattices. Our construction is an extension of Construction 
A and hence we refer to it as Construction A'. We now define polynomial rings and codes over 
polynomial rings. 

Definition 3: (Ch. 4 in [fT6l ) We define the polynomial quotient ring Um = F2[m]/m™ in 
variable u for any m > 1 as 



m— 1 



W™ = < ^ bku'' mod m" | 6^ G F2 

I fc=0 

with regular polynomial addition and multiplication over F2 coefficients along with the quotient 
operation m™ = 0, which is equivalent to cancelling all the terms of degree greater than or equal 
to m. 

Definition 4: A linear code C over Um is a subset of W^ which can be obtained through a 
generator matrix G G W^^" as 

C = {zG I Vz e U^}, 



'Unlike uncoded communication, gray-mapping on Z[i]/piZ[i] is not necessarily optimal since it does not guarantee that the 
neighbouring lattice points in the lattice code are separated by maximum number of information bits. Efficient bit labelling of 
lattice codes is a separate problem of its own and is out of the scope of this work. 
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for some k < n and the matrix multiplication is over the ring Um- 

III. Construction A' of BW Lattice 

We now introduce Construction A' in the following definition. 

Definition 5: A complex lattice A is obtained by Construction A' from a linear code C over 
Um for some m > 1 if A can be written as 

A = <^{u"')Z[i]"' + SC, (5) 

where EC = {$(c) | Vc G C} C Z[i]" is a lattice code obtained from the linear code C through 
the mapping $ : Um -> '^[i] given by 

(m— 1 \ m—1 

j=0 J j=0 

such that ip -.¥2 -^ Z[i] given by ip{0) = and ^^(l) = 1, and $(m) = 1 + i. 

Note that Construction A can be obtained as a special case from Construction A' when m = 1, 
wherein the embedding operation $ coincides with ip given in Definition [2l In the following 
subsections, we use Construction A' to obtain complex BW lattices of dimension 2™ for any 
m > 1 by embedding a linear code C (denoted by C2"0 over the quotient ring Um to a lattice 
code £C (denoted by £C2'n). 

A. Linear codes for construction A' 

In order to obtain BW2^^ as Construction A , we first need to find a suitable linear code €2^ 
over the ring Um- We propose such a linear code which can be obtained by the following the 
generator matrix 



G, 



"1 (sm 
1 1 



u 
where the tensor operation is over the ring Um- 
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Example 1: To obtain 5W4, the linear code C4 can be generated using 
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1 


1 





u 





u 








u 


u 















G U^''^ 



Encoding of linear code €2^ 

By using G2"i as a matrix over Um, the code 62^ is obtained as follows: Let z G W^"", i.e., 
the j-th component of z is given by 



m—l 



h = Yl ^^J^^' 



(6) 



fc=0 



where 6/cj G F2 for all k,j. Using z and G2"i, the code €2^^ C W^'" can be obtained as 

Cam = {x = ZG2- I Vz G W^" } , 



(7) 



where the matrix multiplication is over U^. 

We now provide an example for the proposed encoding technique, showing the positions of 
the information bits that get encoded to the codewords of Cam. 

Example 2: For m = 2, the input vector z and the generator matrix G4 are of the form. 



^0,1 + hi,iu 




1 


1 


1 1 


bo,2 


and G4 = 





u 


u 


bo,3 










u u 

















We define the rate of the linear code Cam as the ratio of the number of information bits per 
codeword and the length of the code (which is also known as the spectral-efficiency of the code). 
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Proposition 1: The rate of the code ^2^ is y. 

Proof: Each component of z carries m information bits in the variables 6/,. j as shown in 
^. This amounts to a total of m2™ bits carried by z. However, since the matrix multiplication 
is over Um, not all the information bits hkj are encoded as codewords of C2m (since n'^ = for 
k > m). Using the structure of G2™ it is possible to identify the indices {k,j) of information 
bits 6fc j which get encoded into the codewords of C2m as follows. Let the set Xg denote the 
indices of the rows of G2m whose components take values or u'^ for each g = 0,l,...,m — 1. 
Due to the quotient operation m™ = 0, the components of z which are in the index set Xq are 
restricted to be of the form, 



m—l—q 



^3 



Y, h.u'yjelq 



fc=0 

For example, zi = YlT=o ^k,iu^ and Z2^ = 0. Using the structure of G2™ we observe that the 
cardinality of Xq denoted by |Xg| is C^, and hence we find the total number of information bits 
per codeword of C2- as EL'o^l^ - k)C]^ = 12"^. ■ 

We now show the equivalence of our encoding technique to Construction D. In other words, 
the following theorem shows that the codewords generated in (|7]) can be uniquely represented 
as vectors of a multi-level code of nested RM codes. 

Theorem 1: The codewords generated in (|7]) can be uniquely represented as codewords ob- 
tained through Construction D. 

Proof: See the proof of Theorem 1 in [|24l. ■ 

Till now, we have presented the linear code €2^ and its encoding technique over the quotient 
ring Um- Now, we discuss the embedding operation of C2m into the Euclidean space. By using 
the map $(n) = 1 + z on C2m, we get the lattice code SC2m. Note that SC2"^ can be used 
as a tile in constructing the BW lattice, i.e., BW2^^^ can be obtained by replicating £62^ in 
Z[i]^™ as BW2m = (1 + z)"'Z[z]^'" + £^C2m. It can be verified that £C2-^ is an arbitrary subset 
of BW2m and does not have cubic shaping. In Fig. [H we plot the complex points generated as 
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Tiling in Barnes-Wall lattices using a lattice Constellation for m = 10 
60 r 




Fig. 1. Filling the complex plane using the tile generated by X^^y (1 + iYbr for m = 10. 



{Er=o (1 + ^y^r I 6,. e {0, 1}} for m = 10. Note that the points generated by ^"Xoi^ + iYbr 
are marked in black, whereas the points in other shades correspond to the shifted version of 
E^~o^(l + ^ybr by constants (1 + z)™, z(l + z)'" and (1 + i)(l + z)™. 

Note that the code £€2^^ does not have good shaping, we observe that the average transmit 
power of the scheme is not small. To fix this problem, we propose a one-to-one mapping (p on 
£C2m to obtain a new lattice code denoted by £2™ such that it has good shaping property. 



B. BW lattice codes with cubic shaping 

Here, we propose a one-to-one mapping on £62^ to obtain a new lattice code £2™ which 
has the cubic shaping property when m is even, and the rectangular shaping property when m 
is odd. For any x = [xi,X2,a;3, . . . ,X2'"] G £€2^^, the mapping operates on each component 
of X as. 



(f){Xj) 



Xj mod 2 2 , when m is even; 
(p (xj mod 2^^ j , when m is odd, 
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where ip{-) is defined on Z m+i [i] as 

z, when 5(z) < 2~2— ; 
z + (2^ - 12^") , when 5R(z) < 2^ 
f{z) = { and 3(z) > 2^; 

z - [2^ + 12^") , when 5R(z) > 2^ 
and 3(z) > 2^. 
Tlie mapping guarantees the following property on £2^' 

|Z2^[i]} , if m is even; 



(9) 



(10) 



r 1^'" r 1^'" 

< Z m+i ^ + 2 < Z m-1 ^ , if m is odd. 
From (flOl) . note that each component of the vector in £2'" is in a cubic box and a rectangular box, 

when m is even and odd, respectively. In Fig. [2l we present the complex points YlT=f) (-'- '^'^Y^r 

with and without the mapping for m = 10. With this, the lattice code £2™ can be obtained 

from C2m through the composition map 

X = 0(<^(-)), (11) 

where $ and (p are given in Definition |5] and ([8]) respectively. The following proposition shows 
that x(-) is a one-to-one map on C2™ 

Proposition 2: The mapping x given in (fTTl) is one-to-one. 

Proof: Since x is a composition mapping of $ and 0, and $(■) is a substitution operation 
using binary representation of complex numbers over the base (1 + i), we have to prove that (p 
given in ([8]) is one-to-one. Here, we provide the proof when m is even. For any Xi,X2 G SC^^ 
such that Xi 7^ X2, we prove that 0(xi) 7^ 0(x2). Applying the modulo operation in ([8]), Xj 
satisfies Xj = 2^rj + 0(xj) for each j = 1, 2, where 0(xj) G £2™ and r^ G 'L[iY'\ This implies 

0(x,) = X, - 2^rj = X, + (1 + z)'"r;, (12) 

for some r^ G Z[z]^™. The second equality follows as (1 + z)™ = a2^, where a G {1, — 1, z, — z}. 
Further, since each component of Xj is of the form Yl^=oi^ + "^YK for 6,. G {0, 1}, the R.H.S 
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lattice constellation with and without cubic shaping 



30 



20 



-20 L 



Fig. 2. Complex points generated by X^^o (^ + iYbr and (j}{J2T=o (^ + iYbr) for m — 10. 



of (fT2l) is nothing but the binary decomposition of (pixj) over the base (1 + i). Since the radix 
representation over (1 + i) is unique, we have 0(xi) = (/)(x2) only if Xi = X2. This completes 
the proof when m is even. The one-to-one nature of can be proved on the similar lines when 
m is odd. ■ 

The above proposition implies that mapping provides a new lattice code with better shaping 
property. The following theorem shows that £2^ can be used as a tile to obtain BW lattices. 

Theorem 2: The lattice code £2™ and the lattice BW2^n are related as BW2^ = (l+z)'"Z[z]^'"© 

L,2'"- 

Proof: See the proof of Theorem 2 in [|24l . ■ 

Using the results of Theorem |2l BW2m is given by BW2m = (1 + 2)'"Z[?]^'" © £2"^, where 
£2™ is the lattice code obtained from C2™ through the mapping x = 0('^(")) ^^ ^m- 



IV. On the Error Performance of the SBWD 

In this section, we study the error performance of the SBWD in decoding the infinite BW 
lattice. In lfT3l . it is shown that for x G BW2m, if there exists y G C^™ such that (i,^j„(x, y) < ^, 
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where N = 2™, then the SBWD correctly finds (or decodes) the lattice point x = x. In the 
context of using SBWD in AWGN channels, the vector y corresponds to y = x + n, where 
X G BW2m and Uj ~ CA/'(0,(t^) Vj. This implies that the codeword error rate (CER) of the 
SBWD given by Pr(x ^ x) is upper bounded as 



Pr(x 7^ x) < Pr ( |np > 

Note that ^ is the packing radius of BW2"^, and hence the above bound is the well known 
sphere upper bound (SUB) [[23 • In |[T3l . the focus was only on the complexity of the decoder but 
not on the analysis of the tightness of the SUB. In other words, the possibility of correct decision 
is not known when |np > ^. We study the error performance and show that the decoder is 
powerful in making correct decisions well beyond the packing radius. Without loss of generality, 
we study the error performance when the zero lattice point is transmitted. We analyze the SBWD 
algorithm and point out the reason for the improvement in the error performance (with reference 
to the SUB). We first recall the SBWD algorithm of [[T3l . 

The Sequential BW Lattice Decoding Algorithm: 



function SEQBW(r, y) 




if y G C^ and A^ < 2'' 




return [yj; 




else 




b = \^iy)\ + r^(y)J mod 2; 




p=l-2(max(|[gf?(y)J-gf?(y), \^{y)\\ 


-s(y))); 


c = RMDEC(r,b,p); 




v = SEQBW(r + l,(y-c)/(l + z)); 




return c + (1 + i)y; 




end if 





January 9, 2013 DRAFT 



16 



end function 



The above decoder is a successive interference cancellation (SIC) type decoder which exploits 
the BW lattice structure as a multi-level code of nested RM codes (as per Construction D). 
At each level, the algorithm uses a variant of the soft-input RM decoder ETI (denoted by the 
function RMDEC which is given as Algorithm 3 in |fT3l ) to decode, and cancel the RM codeword 
at that level. Therefore, the error performance of the SBWD is fundamentally determined by the 
error performance of the underlying soft-input RM decoders . In particular, we have 

Pr(x^x) = Pr(|j£(c, ^c,)j, (13) 

where £{Cr ^ Cr) denotes an error event while decoding TZAi{r,m). Hence, it is important to 
compute Pr(Cr. ^ c,) for each 7lAi{r,m). Along that direction, it is necessary to model the 
effective binary channel induced for each RM code TZAi{r, m). We propose a model for such a 
binary channel which is accurate for r = 0, while for r 7^ 0, it neglects the error propagation in 
the SIC decoder algorithm. To decode the RM code at each level, a hard-decision binary value 
hj is obtained from yj as 

bj = my,)\ + rS(%)J mod 2. (14) 

Due to the combination of the round and the modulo operation (henceforth referred to as the 
round-modulo operation) in (IT4] |. the codewords of TZAi{r,m) are passed through a virtual 
binary channel with the cross-over probability given by, 

P, = Fribj = 1 I c, = 0), 

where c G TZAi{r,m). Since the zero lattice point is transmitted, c is the all zero codeword 
for each TlAi{r,m), and hence the relevant cross-over probability is Pr(6j = 1 \ Cj = 0). The 
following theorem shows that Pc can be upper bounded by a Jacobi-Theta function [|26l . 
Theorem 3: The cross-over probability P^, induced by the round-modulo operation in ([14] ) is 
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Fig. 3. Comparison of the cross-over probability with the upper bound using the Jacobi-Theta function 



upper bounded as 



Pc< e-i^h9 



i4 



9 ' 9 



(15) 



where •§ (z, r) is the Jacobi-Theta function given by 

oo 



TTJa T+27ria2 



Proof: We first compute Pc, and then propose an upper bound. To assist compute Pc, we 
compute the probability that 'St{yj) (or '^{yj)) falls within an interval (z — 0.5, z + 0.5] centred 
around an integer z, when Cj = 0. Since the additive noise is circularly symmetric, it is sufficient 
to calculate the above probability for either ^{yj) or '^{yj). We use y to denote either ^{yj) or 
Q{yj)- For the odd integer case, we have 

oo 

Po = J] Pr(2a + 0.5<i/<2a+1.5), 



a=— oo 

oo 

E 

a=— oo 

oo 

E 



2a+1.5 



^y{y)dy 



Q 



2a+0.5 

2a + 0.5 



a 



/V2 



Q 



2a + 1.5 



a 



/V2 



(16) 
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2 

where Py(y) is the probability density function of y, Q{x) = -^ J^ e~^du, and cr^/2 is the 
variance of y. For the even integer case, we have 

oo 

Pe = ^ P^ (2a - 0.5 < 2/ < 2a + 0.5) , 

f2a+0.5 



a=— oo 

oo 

E 

a=— oo 

oo 

E 



Py{y)dy 



lJ2a-0.5 

, 2a -0.5 



Q 



2a + 0.5 



(17) 



a/V2 J V f^/v^ 

Note that bj is 1 whenever \'Si{yj)\ + \^iyj)\ is an odd number. This can happen when (?) 

l^iUj)! is odd and [Q'(2/j)J is even, or (ii) \^{yj)\ is even and [Q'(2/j)J is odd. From (fT6] ) and 

(fTTl) , we can write 



Po(l-Po) + (l-Po)Po, 



= 2Po-2(Po)'. 
By dropping the term 2(Po)^, we upper bound Pc as 



(18) 



(19) 



Pc < 2Po, 



< 



2E 

a=— oo 

oo 



a=— oo 

(0-5)^ 



Q 



2a + 0.5 



0" 



/v^ 



^3 



E 



(20) 
(21) 



e <^^ 



e i^ ] i7 



9 ' 9 / ' 



where the bound in (l20l) comes from dropping the terms of the form Q ( ^^jj^ ) in (fT6l) . and 



the bound in (|2T1) is due to the Chemoff bound Q{x) < ^e^~. ■ 

Note that the Jacobi-Theta function can be evaluated at any pair (r, z). In Fig. [3l the empirical 
values of Pc are presented along with the bound in (fTST ) for various values of SNR = -^. We 
point out that the bound is not tight due to the Chemoff-bound on each Q(-) function. 
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It is well known that Pc determines the error-performance of a hard decision decoder. Since 
we have a soft-input decoder, we need to obtain the relevant statistics on the soft inputs. We now 
study the soft-input RM decoder used in the SBWD. Unlike the codewords of RM code in ll2T]| . 
the RM codewords at each level of BW lattice take values over {0, 1}. The soft-input used for 
the RM decoder is p = 1 -2d, where d = max (| [3f?(y)J - 3fi(y)|, | [S(y)J - S(y)|). Also, unlike 
the soft metric in [|2T1l . pj is bounded in the interval [0, 1]. This is because dj G [0,0.5], which 
is a result of the round-modulo operation in (fT4] ). One could imagine b and p to be obtained 
from the received vector in a virtual additive noise channel, wherein each component of the 
received vector is always within a distance of 0.5 from either or 1. Therefore, if c denotes a 
RM codeword at a particular level of the transmitted BW lattice point, then the effective noise 
n'^-f'^ as seen by the soft-input RM decoder at that level is of the form, 

dj, when bj = c,; 
nf=l ' ' ' (22) 

I I — dj, when bj ^ Cj] 

for 1 < j < A^. Note that rij has bounded support in the interval [0, 1]. For an analogy with 
respect to the model in llITl . the code alphabet {0, 1} in lfT3l corresponds to the code alphabet 
{-1, 1} in dm and the effective noise n^J'^ in [[HI corresponds to the AWGN in [|2l]- At each 
level of the BW lattice, the lattice code (1 + iyTZAi{r,m) for any < r < ??7 — 1 has the 
minimum squared Euclidean distance of N. By using the proposition in Section IV. A of [2T|, 
the probability of incorrect decision of the soft-input RM decoder at each level of SBWD is 
upper bounded as shown in the proposition below. 

Proposition 3: The codeword error rate Pr(Cr 7^ c^) for each TZM.{r, m) is upper bounded as, 

Pr(c^ ^c,.) <Pr ( |n"-^^P > — J for r = 0, 1, . . . ,m - 1. (23) 

It is important to note that the above bound is different from Pr(|iip > ^) since n is Gaussian 



distributed. We do not have closed form expression on the distribution of either rij or 



n. 



In Fig. HI we display the histogram of the realizations of n'^-^-^ for various values of cr^, when 
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Fig. 4. Histogram of nf^ when SNR = ^ takes the values dB, 5 dB, 10 dB and 25 dB. 



the zero RM codeword is the transmitted. Note that for cr^ = dB, the histogram of Uj has 
the triangular shape centred around 0.5, which implies a very high (close to 0.5) cross-over 
probability when obtaining the hard decision vector b. On the other hand, at lower values of a^, 
the distribution is skewed towards zero indicating smaller cross-over probability. 



V. SBWD TO Decode BW Lattice Code £om for AWGN Channel 



In this section, we discuss the use of SBWD to decode the lattice code £2™- First, we describe 
a method to transmit the codewords of £2™- For any x G £2™, the transmitted vector is of the 
forme 



xt = (2x - c) 



(24) 



The transmitted vector is offset by a constant c towards the origin to reduce the average transmit energy. 
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where 

(2^ — I) + i (2^ — 1) , when m is even: 
( 2^^ — l) + i ( 2^^ "-'-)' when m is odd. 
Using the scale and the shift operation in (l24l ). each component of x^ takes value from the regular 
2"*-QAM constellation. In particular, the QAM constellation is square and non-square when m 
is even and odd, respectively. When x^ is transmitted, the received vector y is given by 

y = Xi + n, (26) 

where n is the AWGN with Uj ~ CA/'(0, a^) Vj. In this section, SNR of the channel is defined 
as Eg/a'^, where Es denotes the average energy of 2" -QAM constellation. With the inverse 
operation to (|24l) as y = ^y + c, the equivalent AWGN channel becomes 

y = X + n, (27) 



^2^ 



where x G £2- and Uj ~ CA/'(0, ^). We use the SBWD O on ^ to decode the lattice code 
£2'"- When a codeword of £2™ is transmitted, the SBWD decodes to a lattice point in the infinite 
lattice BW2rr^ ■ In such a decoding method, irrespective of whether the decoded lattice point falls 
in the code or not, the information bits can be recovered from the decoded RM codewords at 
every level of SBWD (as shown in the algorithm in Sec. HVT ). 

A. Simulation results on the codeword error rate (CER) of SBWD 

In this subsection, we present the CER of the SBWD along with some upper bounds and lower 
bounds. For the simulation results, we use SNR = E^/a'^, where Eg denotes the average energy 
of the regular 2™-QAM constellation. In each of Fig. [SHU we present (/) the CER of the SBWD, 
(//) the SUB (Section IV.D, [22l), (Hi) the sphere lower bound (SLB) (Section IV.D, [|22l), (iv) 
the CER in decoding TZAi{0,m) at the first level of the SBWD, and (v) the upper bound on 
the CER in decoding 7^A^(0, m) given by Prdn'^^-^p > ^) (obtained through simulation results 
by empirically generating n"^-^-^). 
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Fig. 5. CER of SBWD decoding BW4. 
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Fig. 6. CER of SBWD for decoding BWw. 



From Fig. |5]|9l we make the following observations: the SUB is not a tight upper bound on 
the CER of SBWD. Also, Pr(|n^^-^|2 > f ) is an upper bound on the CER of SBWD, and in 
particular, it is a tighter upper bound than the SUB. The CER of the soft-input RM decoder 
for 7^A^(0, m) is a tight lower bound on the CER of the SBWD. This implies that if there is 
no error at the first level of the decoder, then with high probability, there will be no errors at 
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Fig. 7. CER of SBWD for decoding 5^64. 
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Fig. 8. CER of SBWD for decoding BW256- 

subsequent levels of the soft-input RM decoder. In summary, the simulation results highlight that 
the SBWD is quite powerful in making correct decisions even beyond the packing radius, and 
the deviation from the SUB increases for larger dimensions. As a result SBWD can be employed 
to efficiently decode lattice codes of large block lengths with low-complexity. This behaviour in 
the error performance of SBWD was not known in the literature. 
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Fig. 9. CER of SBWD for decoding BWio24. 



B. Comparing the complexity of the SBWD with the list decoder in /[75]/ 

In this subsection, we compare the complexity of the SBWD with the BW list decoder lITSl . 
For a fair comparison, we assume that the list decoder is implemented on a single processor. 
On a single processor, the complexity of the SBWD is 0(Mog^(A^)), whereas the complexity 
of the list decoder is 0(A^^)(/(m, r/))^, where l{m,r]) is the worst case list size at a relative 
squared distance of r] (the relative squared distance is the squared Euclidean distance normalized 
by the dimension of the lattice). We compare the complexity of the two decoders for a codeword 
error rate of 10^'^. In particular, we first approximate the error performance of the SBWD as a 
bounded distance decoder for some radius f/, and then compute the complexity of the list decoder 
with the corresponding value of fj. In Table HI we display the lower bound (as given in Theorem 
1.3 in [[TSl ) on the complexity of the list decoder to achieve the error performance of SBWD. 
The table shows that the list decoder has higher complexity than the SBWD to achieve the same 
performance. In summary, for single processor implementation, SBWD can be preferred to the 
list decoder to decode BW lattice codes of large block lengths. However, for codeword error 
rates lower than that of SBWD, the list decoder has to be used, preferably on parallel processors. 
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TABLE I 
Complexity of the list decoder ifTSll to achieve the performance of SBWD 



Dimension N 


V 


A lower bound on N^{l{m,7]))'^ 


Mog^(Ar) (complexity of SBWD) 


4 


0.33 


16 


16 


16 


0.4 


256 


256 


64 


0.48 


4096 


2304 


256 


0.56 


262144 


16384 


1024 


0.67 


1.07 X 10^ 


102400 



Table H also shows the potential of SBWD to decode well beyond the relative squared distance 
of ?7 = 0.25. For complex dimensions of 256 and 1024, the effective radius of SBWD is as high 
as y and ^, respectively. 



VI. Noise Trimming Technique for the SBWD 

When a codeword of £2™ is transmitted, the SBWD decodes to a lattice point in the infinite 
lattice BW2m. In such a decoding method, irrespective of whether the decoded lattice point 
falls in the code or not, the information bits can be recovered from the decoded RM codewords 
at every level of SBWD (as shown in the algorithm in Sec. HVl) . To further improve the error 
performance, we force the SBWD to specifically decode to a codeword in the lattice code, and 
subsequently recover the information bits, with more reliability. We refer to such a decoder as 
the BW lattice code decoder (BWCD). We use a technique that forces the SBWD to decode to 
a codeword in the lattice code £2™ • We refer to this technique as the noise trimming technique, 
which exploits the structure of £2™- From (fTOl ). we know that each component of a codeword 
is within a rectangular box i3 C C. In particular, the box B shares its edges with Z^^ [i] and 
Z m+i + zZ m-i when m is even and odd, respectively. In order to use SBWD, and to decode 
to a codeword within the code, we trim the in-phase and quadrature components of the received 
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vector (the algorithm is given below) to be within a box B' ^ B marginally larger than B by 
length e on each dimension. Then, we feed the trimmed received vector to the SBWD and decode 
the information bits. Note that the choice of e is crucial to decode a codeword within the code, 
and to improve the BER with reference to the SBWD. We now provide an algorithm for the 
trimming method, which works independently on the in-phase and quadrature component of the 
scalars in y = [yi,y2, ■ ■ ■ , 1/2™] in (l27] ). In particular, the algorithm presented in the sequel works 
on the in-phase and quadrature component of yj when m is even. Extension to the case when 
m is odd is straightforward. 

Algorithm for the trimming technique when m is even: 



Input y eR (either ^{yj) or ^{yj)) 
function TRIM(y, e) 
A = (2f -l)/2 
r = y - A 
t = A + e 
if \r\ > t 
s = t/\r\ 
b = s X r 
else 

b = r 
end if 

return b + A 
end function 



Using BWCD, we have obtained BER for dimensions when m = 2,4, and 6, and compared 
them with the BER of the SBWD. The plots as shown in Fig. [TOl indicate that BWCD outperforms 
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SBWD by 0.5 dB. For the presented results, we have used e = -^, which corresponds to the 
packing radius of ^. The above value of e was optimized based on the simulation results 
by comparing the BER for various values of e. Intuitively, trimming the received vector to fall 
within the packing radius of a lattice point in the edge of the lattice code forces to SBWD to 
decode to a lattice point in the edge of the code rather than a lattice point outside the lattice 
code. 




Fig. 10. Comparison of BER between BWCD and SBWD for m = 2, 4, and 6. 



VII. Conclusion and Directions for Future Work 

In the first part of this paper, we have introduced a new method of encoding complex BW 
lattices, which facilitates bit labelling of BW lattice points. As a generalization, the proposed 
technique is applicable to encode all Construction D complex lattices. In the second part of 
this paper, we have used complex BW lattice codes for communication over AWGN channels. 
To encode the code, we have used Construction A', and to decode the code we have used the 
SBWD. We have studied the error performance of the SBWD, and have shown that the Jacobi- 
Theta functions can characterize the virtual binary channels that arise in the decoding process. 
We have also shown that the SBWD is powerful in making correct decisions beyond the packing 
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radius. Subsequently, we have used the SB WD to decode the complex lattice code through the 
noise trimming technique. This is the first work that uncovers the potential of SBWD (in terms 
of the error performance) in decoding lattice codes of large -block lengths with low-complexity. 
This work can be extended in one of the following ways: 

• The SBWD proposed in [fT3l uses a soft-input, hard-output RM decoder at each level of 
Construction D. It will be interesting to study the error performance of the lattice decoder 
with soft-input, soft-output iterative RM decoders. 

• We have presented the error performance of the SBWD through simulation results, and 
hence we now know the SBWD error performance with reference to the sphere lower 
bound and the sphere upper bound. A closed form expression on the error performance of 
the SBWD could be obtained for a better understanding of the decoder performance. 
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